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We numerically calculate two-particle correlation functions of CERN-SPS 158 A GeV 
Pb+Pb central collisions based on a (3+l)-dimensional relativistic hydrodynamics with 
a first order phase transition. We analyze the pair momentum dependence of pion source 
sizes extracted from the Yano-Koonin-Podgoretskh (YKP) parametrization which is ex- 
pected to give the source sizes directly We find that, even in the case of the first order 
phase transition, the collective expansion and surface dominant freeze-out of the fluid 
naturally lead to the opaque source for which the interpretation of the temporal source 
parameter as the emission duration breaks down. 

1. INTRODUCTION 

In the physics of relativistic heavy ion collisions aiming at studying the QCD phase 
transition of hot and dense matter, two-particle correlation is widely used to investigate 
the characteristics of the space-time evolution in the collision processes [[!]]. It is well 
known as Hanbury-Brown Twiss effect that the two-particle intensity correlation function 
contains the information on the "size" of the source of emitted particles due to a quantum 
symmetry. In recent experiments, the sizes are estimated from the correlation function 
through Gaussian fitting and the physical meaning of the sizes depends on the fitting 
function. For an azimuthally symmetric source, the Yano-Koonin-Podgoretskii(YKP) 
parametrization is used as one of the fitting functions. In this parametrization, three 
of fitting parameters can be regarded as source sizes directly. This advantage gives us a 
chance to observe a signature of quark-gluon matter because large emission duration is 
expected in the case of a first order phase transition In the case of the relativistic heavy 
ion collisions, the reaction is highly dynamic; the particle source cannot be considered to 
be static one. So we analyze the size parameters for pions quantitatively based on a 
hydrodynamic model which takes into account the longitudinal and transverse expansion 
and a first order phase transition ||. 
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2. FORMALISM 

From a hydrodynamic point of view, we consider the pion source as a thermalized one 
and a completely chaotic source. Hence, the two-particle correlation function is given as 

C(k 1; k 2 ) = 1 + {dN }^yx /( p ka y (!) 

where ki and k 2 are momenta of a emitted pion pair and J(ki,k2) is the interference 
term. 

As usual, we introduce the relative momentum q^ = kf^ — kf^ and the average momentum 
= (ki + k%)/2 which can be put as = (K°, K?, 0, Kl) by virtue of the cylindrical 
symmetry. Due to the on-shell condition, q^K^ = 0, only three components of the relative 
momenta are independent. In the YKP parametrization, q±_ = Jq\ + q 2 , q\\ = q z , and q° 
are used as the independent components. Then, the fitting function is given by 

C(q^K») = 1 + Aexp {-R\{K^)q\ - R\{K»)[<§ - (g ) 2 ] 

-[i2 2 (^)+i?J(^)][ ? X(^)] 2 }> (2) 

where w M is a four- velocity which has only longitudinal component. The size parameters 
to be determined by a fit are RiS (i = _L, ||, 0). We may choose a special reference frame 
where the velocity parameter v vanishes (YKP frame), because the three size parameters 
are invariant under longitudinal boosts. Introducing the source function, 

and the weighted average, 

_ fMg)W 
{A{X }) " Jd*xS(x»,K») ' (4) 

the size parameters in this frame (YKP frame) are expressed as 

Rl(K») = (Ayf, (5) 
RliK^ = (Atf - ^(xt) + ±[(Axf - (Ay)% (6) 

PL P± 



f, {xz)+ i 



R\{K») = (Az) 2 -^{xz) + ^[(Ax) 2 -(Ay)% (7) 



where x = x - (x), Ax = ^ {x 2 ) - (x) 2 , (3 ± = K T /K° and /3\\ = K L /K°. 

In order to regard the above R^s as naive source sizes directly, the second term and the 
third term in the right hand side of Eqs. (@) and (0) must be small compared to the first 
term of each equation. Though it has been shown that this condition holds within a class 
of thermal model @, it is not trivial in other cases. The case of a hydrodynamic model 
as a more realistic one is the point of the present paper. In the following, we calculate 
three quantities: HBT radii extracted from the fit to Eq. ([]]), the space-time extensions, 
Eqs.(|)-(0), and the source sizes, the first terms in Eqs.® and (0). 
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3. DISCUSSION 
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Fig. 1 Kt dependence of YKP radii. Open and closed circles show the ex- 
perimental data and HBT radii, respectively. Solid lines stand for space-time 
extensions (||)-(0). Dotted lines stand for source sizes. (Az for R» and At for 
Ro.) 

Figure 1 shows Kt dependence of the pion HBT radii. Close circles stand for the HBT 
radii. Solid lines are space-time extensions which are expected to agree with the HBT 
radii when only thermal pions are considered. Dotted lines stand for the source sizes (Az 
for and At for R ). As far as R± and Ru are concerned, those quantities seem to be 
consistent with one another and well reproduce the NA49 experimental results Q5| (open 
circles). On the other hand, Rq is smaller than the experimental one in spite of the fact 
our calculation includes a first order phase transition. The time durations At agrees with 
Rq at small Kt- However, some deviations from R can be seen at large Kt region, that 
means a failure of interpretation of Rq as the time duration. We find that this difference 
occurs from the third term in Eq.(^). 

In Eq.(HD, the two additional terms means correlation between x and z (the second 
term) and the thickness of the source (the third term). Though contribution from the 
second term is small in our model, the third term becomes negative in same order with 
(At) 2 ||. Such a kind of the source is called an "opaque source". The "opaque" source 
emits particles dominantly from the surface and the "transparent" source emits from 
whole region. Figure 2 shows the source function @ projected onto the x — y plane, 
St(x, y) = * J dzdtS(x fJl , K^) being normalized as J dxdyST(x, y) = 1. Note that the 
average transverse momentum of the emitted pions is parallel to x direction in the figure. 
The left figure shows clearly surface emission; pions are emitted from the crescent region 
and the source is thinner in the x direction than in the y direction. This is the typical 
property of the opaque source which appears in the factor (Ax) 2 — (Ay) 2 in Eq.(^). Note 
that this factor also appears in the longitudinal components of the size, Eq.(^), but the 
prefactor (3 2 makes the contribution to R 2 small. If we neglected the transverse flow by 
artificially putting U r = 0, the source function becomes as shown in the right figure in 
Fig. 2. In this case, the source function is proportional to the space-time volume of freeze- 
out hypersurface. As a consequence of restoration of the azimuthal symmetry which is 
shown clearly in the figure, the above measure of source opacity vanishes. Although the 
particle emission takes place almost from the thin surface, the source is not opaque in this 
sense. When the transverse flow exists, the source function is deformed by thermal Boltz- 
mann factor exp(i^T sinh Yt cos <p/Tf) with Yt almost proportional to r. Consequently, 
the number of emitted particles increases in the region where x > (i.e., cos0 > 0) and 
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decreases in the region where x < (cos0 < 0). This flow effect deforms the surface 
dominant volume distribution (right figure) to the crescent shape (left figure). 
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Fig. 2 The source functions as functions of transverse coordinates x and y for 
K T = 450 MeV and Y n7T = 4.15 (Y n7T = \ In gSzf^ )- Th e left figure denotes 
the source function with the transverse flow. The right figure denotes the 
source function without transverse flow. 

In summary, we analyze dependence of source parameters of the YKP parametriza- 
tion based on the relativistic hydrodynamics for the CERN-SPS 158 GeV/A Pb+Pb 
collisions. We obtain the results which are almost consistent with the experiment. How- 
ever, the temporal source parameter, i? , shows smaller value than the experiment. The 
source opacity makes the interpretation of Rq as the time duration doubtful at large K T 
region. We find that the source opacity is caused by the transverse flow and the charac- 
teristics of the freeze-out hypersurface, the surface dominant freeze-out. The deviations 
of our results from the experiment would be improved by including the resonance decay 
and other effects. 
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